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Figure  V,  A  schematic  representation  of  the  propagation  processes  which  need  to  be 
included  in  the  description  of  mantle  phases  (after  Kennett  1989b). 

Figure  2;  Script  file  display  for  the  adaptation  of  the  iasp91  travel  time  routines  to 
display  phase  and  group  velocities  in  order  to  compare  the  ray  based 
assignments  of  conventional  location  procedures  with  the  discrimination 
phases  Pg,  Lg. 

Figure  3:  WKBJ  synthetic  seismogram  for  a  simple  crustal  model  representative  of  a 
shield  region.  Up  to  4  crustal  multiples  are  allowed  for  reflection  from  each 
interface  or  for  diving  waves.  The  P  and  S  wave  calculations  were  carried  out 
separately  and  superimposed  in  plotting. 

Figure  4:  Recordings  of  an  earthquake  (FF)  and  a  quarry  blast  (SM)  along  the  same 
array  of  portable  seismographs  in  "southeastern  Australia.  In  each  case  the  Pn 
phase  can  be  clearly  observed  even  though  the  earthquake  phase  are  dose  to  a 
node  for  P,  it  is  much  more  difficult  to  detect  the  Sn  arrival  even  with  a 
number  of  records. 

Figure  5:  P  and  S  wave  sections  for  the  vertical  component  of  displacement  for  an 
earthquake  in  northwestern  Scotland  recorded  during  the  long-range  refraction 
profile  through  the  British  Isles.  The  reduction  velocity  for  the  P  section  is  8 
km/s,  and  for  the  S  section  it  is  4.62  km/s  (after  Kaminiski  et  al  1974). 

Figure  6:  Three  component  record  sections  for  a  large  explosive  shot  off  the  French 
coast  recorded  across  a  long  range  profile  in  France  (after  Hirn  1977).  Reduction 
velocity  for  the  Z  section  is  8  km/s,  whereas  for  the  L  and  T  sections  the 
velodty  is  4.62  km/s 

Figure  7:  Ray  diagram  illustrating  the  way  in  which  a  mantle  structure  with  a 
number  of  gradient  zones  can  give  rise  to  a  sequence  of  phases  returned  from 
the  mantle  with  an  apparent  en-enchelon  behaviour  in  t«me. 

Figure  8:  Configuration  of  the  SVEKA  and  BALTIC  and  POLAR  refraction  experiments 
in  Fennoscandia. 

Figure  9:  Comparison  of  P  and  S  record  sections  from  the  SVEKA  refraction  profile  in 
Finland  using  vertical  component  seismometers,  the  traces  are  amplitude 
normalised.  The  reduction  velocity  for  the  P  section  is  8  km/s,  and  for  the 
S  section  it  is  4.62  km/s  (after  Luosto  et  al  1984,  Grad&  Luosto  1987). 

Figure  10:  Comparison  of  P  and  S  record  sections  from  the  BALTIC  refraction  profile 
in  Finland  using  vertical  component  seismometers,  the  traces  are  amplitude 
normalised.  The  reduction  velocity  for  the  P  section  is  8  km/s,  and  for  the 
S  section  it  is  4.62  km/s 
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Figure  11:  Comparison  of  P  and  S  record  sections  from  the  POLAR  refraction  profile 
in  northern  Fennoscandia  using  vertical  component  seismometers,  the  traces 
are  amplitude  normalised.  The  reduction  velocity  for  the  P  section  is  8  km/s, 
and  for  the  S  section  4.62  km/s  (after  Luosto  et  al  1989). 

Figure  12:  Portion  of  a  record  section  assembled  from  hand-digitised  seismograms 
from  a  detailed  DSS  profile  across  the  southern  Ukrainian  shield  with 
reduction  velocity  6  km/s.  The  amplitude  variability  across  the  very  close  trace 
spacing  (around  200  m)  illustrates  the  problems  involved  in  trying  to  use  the 
ratios  of  different  seismic  phases  as  measures  of  source  radiation  patterns  (after 
Jentsch  1979). 

Figure  13:  Portion  of  the  record  section  from  a  DSS  profile  in  northern  Siberia,  the 
reduction  velocity  is  8.2  km/s  (after  Yegorkin  &  Chernyshov  1983). 

Figure  14:  Amplitude  distance  behaviour  for  topography  at  the  crust-mantle 
boundary.  Solid  circles  indicate  calculations  for  the  reference  model  at  the  right 
with  a  Moho  depth  of  32.5  km.  Open  squares  indicate  propagation  from  north 
to  south  and  open  triangles  propagation  from  south  to  north. 

Figure  15:  Record  sections  of  synthetic  seismograms  illustrating  the  effect  of 
topography  at  the  crust-mantle  boundary.  R  reference  structure;  NS 
propagation  from  north  to  south;  SN  propagation  from  south  to  north.  No 
scaling  is  applied  with  range. 


INTRODUCTION 


A  number  of  the  techniques  which  have  been  proposed  for  discriminating 
between  different  types  of  seismic  sources  at  regional  distances  depend  on  the 
relative  amplitudes  of  different  P  and  S  phases.  Although  the  phase  denoted  in 
discrimination  work  as  Lg  is  frequently  the  largest  phase  on  a  regional 
seismogram,  it  is  known  to  be  sensitive  to  variations  in  crustal  structure.  As  a 
result  attention  has  been  transferred  to  Sn  as  a  representative  of  the  S  wave 
radiation  from  the  source,  and  comparison  has  been  made  with  the  Pn  phase  to 
characterise  the  P  radiation. 

However,  detailed  observations  of  the  amplitude  behaviour  of  Pn  and  Sn 
frequently  show  a  very  different  dependence  in  the  range  200  -  300  km  away  from 
the  source.  Whilst  Pn  is  clear  and  separates  distinctly  from  the  rest  of  the  P 
wavefield,  Sn  is  often  not  discernable  until  distances  around  300  km.  Such 
behaviour  has  been  observed  on  refraction  surveys  in  Europe  and  from  regional 
events  in  southeastern  Australia. 

In  order  to  look  at  the  causes  of  this  behaviour  and  the  way  in  which  the  Pn 
and  Sn  amplitudes  vary  with  distance  we  have  examined  a  wide  range  of  long- 
range  refraction  profiles  from  the  Eurasian  area.  We  have  also  investigated  the 
influence  of  topography  at  the  crust-mantle  boundary  on  the  amplitude 
distribution  to  be  expected  from  mantle  phases  such  as  Pn,  Sn. 

THE  PHYSICAL  CHARACTER  OF  REGIONAL  PHASES 

The  analysis  techniques  used  in  discrimination  work  tend  to  label  prominent 
arrivals  on  regional  seismograms  with  simple  descriptors  (e.g.  Pn,  Pg,  Sn,  Lg) 
based  in  large  part  of  the  properties  observed  at  regional  arrays.  However,  in 
order  to  understand  the  way  in  which  the  relative  amplitudes  of  such  arrivals 
change  with  range  we  need  to  work  with  a  more  direct  physical  description  of  the 
propagation  processes  which  give  rise  to  the  classes  of  regional  arrivals. 

Kennett  (1989b)  has  shown  how  a  description  of  the  propagation  processes 
involved  in  the  major  regional  phases  can  be  made  in  terms  of  reflection  and 
transmission  operators  for  propagation  in  a  laterally  heterogeneous  medium. 
Figure  1  sketches  the  classes  of  ingredients  which  have  to  be  taken  into 
consideration  in  the  representation  of  the  mantle  phases  -  both  the  propagation 
characteristics  of  the  crust  and  the  mantle  have  to  be  known  for  effective 
modelling.  In  addition  the  development  of  suitable  computational 
approximations  for  such  operators  is  dependent  on  a  knowledge  of  the  fine 
structure  of  the  seismic  arrivals. 

The  phase  and  group  velocity  definitions  used  for  phases  such  as  Pg,  Lg  do  not 
fit  in  easily  with  the  phase  assignment  schemes  used  by  many  location  schemes, 
which  are  based  on  direct  ray  modelling  of  the  possible  propagation  paths. 
However,  with  only  a  modest  amount  of  effort,  it  has  proved  possible  to  modify 
the  output  of  the  new  iasp91  travel  time  routines  which  implement  the  travel 
time  tables  of  Kennett  &  Engdahl  (1991),  to  provide  specific  values  of  the  phase 
and  group  velocities  associated  with  the  conventional  phase  names.  An  example 
of  a  script  file  generated  by  this  process  is  displayed  in  figure  2.  As  has  been 
previously  noted  by  a  number  of  authors  (see  e.g.  Kennett  1985)  the  definition  of 
Lg  as  a  phase  with  a  group  velocity  of  3.5  km/s  is  appropriate  because  of  the 
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» 1:  A  schematic  representation  of  the  propagation  processes  whi< 
included  in  the  description  of  mantle  phases  (after  Kennett  1989b). 
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Figure  2:  Script  file  display  for  the  adaptation  of  the  iasp91  travel  time  routines  to 
display  phase  and  group  velocities  in  order  to  compare  the  ray  based 
assignments  of  conventional  location  procedures  with  the  discrimination 
phases  Pg,  l.g. 
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mutual  interference  of  many  multiple  S  reflections  within  the  crust,  but  this 
means  that  a  phase  association  via  ray  definition  will  shift  quite  rapidly  between 
different  phases  with  distance.  A  similar  problem  tends  to  happen  with  Sn  as  the 
distance  increases,  the  amplitude  of  the  waves  returned  from  the  mantle 
diminish  and  fall  below  the  coda  of  P,  but  the  surface  multiples  remain  visible 
and  can  be  confused  with  the  true  arrival  Sn  (Kennett  1985).  The  reinforcement 
of  the  surface  multiples  occurs  because  there  are  a  number  of  different  physical 
propagation  paths  which  have  the  same  time  behaviour. 

In  addition,  the  generic  descriptions  Pn,  Sn  used  for  mantle  phases  are  not 
necessarily  an  accurate  description  of  the  character  of  the  arrivals  which  are 
actually  picked  from  a  regional  seismogram.  As  we  shall  see  below,  long-range 
refraction  experiments  ascribe  significant  substructure  to  the  "Pn"  arrival 
between  200  and  800  km.  The  interaction  of  these  multiple  arrivals  from  the 
mantle  and  surface  multiples  gives  a  complex  structure  of  propagation  processes 
which  are  reflected  in  the  visual  complexity  of  regional  seismograms. 

An  inexpensive  and  convenient  tool  for  assessing  the  way  in  which  different 
types  of  propagation  effects  can  interact  to  build  up  the  patterns  seen  on  an 
individual  seismogram  is  provided  by  WKBJ  seismograms  (Chapman  1978)  for  a 
horizontally  stratified  model.  Unlike  methods  based  on  the  full  reflection 
response  of  the  medium,  the  different  classes  of  propagation  paths  have  to  be 
specified  separately.  As  a  result  there  is  a  significant  initial  effort  in  assembling 
the  ray  code  descriptions  which  control  the  class  of  physical  processes  which  are 
included,  but  once  this  is  done  the  computation  of  seismograms  is  very  rapid. 
Further  there  is  no  difficulty  in  isolating  the  contribution  of  different  processes  by 
modifying  the  choice  of  ray  codes.  An  example  of  the  WKBJ  synthetics  for  a 
simple  crustal  model  representative  of  a  shield  region  is  presented  in  figure  3, 
together  with  the  group  velocity  trajectories  associated  with  the  array  definitions 
of  Pg,  Lg.  As  expected  these  trajectories  cut  across  a  sequence  of  multiple  arrivals 
representing  P  and  S  wave  reflections  within  the  crustal  waveguide. 

For  this  WKBJ  example  the  ray  codes  were  chosen  so  that  up  to  4  surface 
multiples  of  each  major  ray  path  would  be  included.  This  includes  reflection 
from  each  of  the  boundaries  in  the  model  or  refracted  (diving)  waves  below 
interfaces.  The  theoretical  seismograms  are  presented  with  trace  normalisation, 
and  clearly  display  the  significance  of  multiples  of  the  mantle  phases  at  greater 
distance.  The  multiple  reflected  waves  in  the  crust  tend  to  asymptote  to  the 
expected  trajectory  of  Lg  and  we  can  see  that  we  would  expect  systematic  errors 
from  the  group  velocity  definition  of  Lg  for  distance  ranges  close  to  the  critical 
point  of  successive  multiples  of  the  reflection  from  the  crust  mantle  boundary. 
As  more  and  more  multiples  contribute  to  the  wavefield  the  association  of  the 
maximum  S  amplitude  with  the  specific  group  velocity  of  Lg  improves. 
However  for  epicentral  distances  less  than  400  km  it  would  be  preferable  to 
attempt  to  associate  the  observed  Lg  arrival  with  the  theoretical  times  for  Sg  or 
SgSg. 

It  may  be  possible  to  introduce  a  partial  allowance  for  lateral  heterogeneity  into 
the  WKBJ  calculation  by  introducing  systematic  perturbations  to  the  central 
phase  functions  which  are  currently  constructed  for  horizontally  stratified  model. 
Such  a  procedure  will  be  investigated  over  the  next  year. 
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Figure  3:  WKBJ  synthetic  seismogram  for  a  simple  crustal  model  representative  of  a 
shield  region.  Up  to  4  crustal  multiples  are  allowed  for  reflection  from  each 
interface  or  for  diving  waves.  The  P  and  S  wave  calculations  were  carried  out 
separately  and  superimposed  in  plotting. 


Observations  of  mantle  phases 


A  number  of  authors  have  noticed  that  the  onset  of  the  Sn  phase  around 
200  km  can  be  difficult  to  pick  by  comparison  with  the  equivalent  Pn  phase  even 
for  those  situations  where  the  P  coda  is  low.  These  observations  come  from  a 
wide  range  of  stable  tectonic  environments.  For  example,  the  phase  association 
procedure  used  for  locating  earthquakes  in  southeastern  Australia  only  declares 
an  Sn  phase  when  the  estimated  distance  is  greater  than  300  km  because  of  the 
difficulty  in  obtaining  reliable  picks  at  shorter  distances.  The  amplitude  of  the 
energy  return  from  the  mantle  is  dictated  by  the  size  of  the  velocity  gradients  for 
P  and  S  waves  in  the  uppermost  part  of  the  mantle.  The  pattern  of  observations 
for  Sn  with  a  low  amplitude  arrival  around  200  km  but  greater  amplitude  at 
ranges  beyond  300  km,  suggests  that  the  S  wave  gradient  is  low  in  the  uppermost 
mantle  and  then  increases  at  a  depth  of  10  or  so  kilometres  below  the  crust- 
mantle  boundary.  Recently  Gajewski  et  al  (1990)  have  shown  that  this  pattern  of 
slow  emergence  of  Sn  is  observed  on  many  refraction  profiles  in  western  Europe 
in  zones  with  a  similar  geologic  setting  to  Eastern  Australia  i.e.  regions  affected  by 
mountain-building  episodes  in  the  last  300  Ma.  They  have  associated  the 
differences  in  velocity  gradient  for  P  and  S  with  a  possible  differences  in  mineral 
composition  with  depth.  Such  geologic  environments  also  occur  in  much  of  the 
Soviet  Union.  As  we  shall  see  later  in  some  areas  of  precambrian  rocks  Pn  and 
Sn  are  observed  with  comparable  clarity,  so  that  there  may  well  be  differences  in 
the  nature  of  the  uppermost  mantle  which  depend  on  the  geological  history  of 
the  region  being  studied.  Gajewski  et  al  (1990)  also  noted  differences  in  the 
character  of  Sg  phases  compared  with  the  corresponding  Pg  arrivals. 

There  are  only  a  limited  number  of  experiments  in  which  both  P  and  S 
propagation  can  be  tracked  to  significant  distances  with  a  station  spacing  close 
enough  for  the  evolution  of  the  wavefield  to  be  judged.  Most  of  these  cases  are 
for  explosions,  but  there  a  very  limited  number  of  earthquake  observations. 
During  the  operation  of  a  linear  array  of  seismographs  in  southeastern 
Australian  during  the  southern  summer  of  1989/1990  we  were  fortunate  to 
record  both  an  earthquake  (FF)  and  a  quarry  blast  (SM)  which  lie  along  the  line  of 
the  profile.  The  corresponding  record  sections  are  displayed  in  figure  4  together 
with  an  inset  map  of  the  array  configuration.  The  general  character  of  the 
wavefield  is  very  similar  although  there  are  differences  in  the  P  radiation  (by 
accident  a  nodal  plane  for  P  wave  radiation  from  the  Fitzroy  Falls  earthquake 
passes  through  the  line  of  the  array).  On  each  section  we  can  see  the  change  in 
phase  velocity  with  the  emergence  of  the  Pn  arrivals  but  the  corresponding 
emergence  of  the  Sn  phase  cannot  be  reliably  tracked.  The  amplitude  of  the 
P  coda  is  quite  high  and  there  is  little  difference  in  frequency  content  between  the 
P  and  S  wave  fields.  This  example  gives  a  good  illustration  of  the  merits  of 
multi-station  recording.  Where  close  station  spacing  was  achieved  there  is  a 
good  correlation  of  the  records  at  successive  and  phases  can  be  readily  identified. 
Even  for  the  isolated  distant  station  identification  of  arrivals  is  simplified  by 
extraction  from  closer  distances.  In  the  first  period  after  the  installation  of  a  new 
regional  station  it  would  be  highly  desirable  to  build  up  a  record  section  from 
available  sources  so  that  the  structure  in  the  neighbourhood  of  that  station  can  be 
calibrated. 
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Figure  4:  Recordings  of  an  earthquake  (FF)  and  a  quarry  blast  (SM)  along  the  same 
array  of  portable  seismographs  in  southeastern  Australia.  In  each  case  the  Pn 
phase  can  be  clearly  observed  even  though  the  earthquake  phase  are  close  to  a 
node  for  P,  it  is  much  more  difficult  to  detect  the  Sn  arrival  even  with  a 
number  of  records. 


A  further  example  of  a  detailed  set  of  earthquake  observations  at  regional 
ranges  is  provided  by  the  recording  of  an  earthquake  in  Scotland  during  the  long- 
range  refraction  experiment  conducted  in  1974.  The  earthquake  was  displaced 
from  the  line  of  the  profile  but  there  is  still  an  very  useful  set  of  observations 
from  70  to  290  km  away  from  the  event  (figure  5).  There  is  a  clear  mantle  S 
arrival  with  a  rather  emergent  onset  which  in  this  case  can  be  tracked  back  to  the 
cross  over  with  the  Sg  phase.  Figure  5  displays  the  onset  of  both  the  P  and  S  wave 
portions  of  the  wavefield  as  recorded  on  vertical  component  seismometers.  The 
travel  times  superimposed  on  the  S  wave  section  are  derived  from  a  model  of 
the  crustal  P  wave  velocity  determined  by  refraction  work  by  assuming  a  ration  of 
P  and  S  velocities  of  V3.  These  theoretical  times  are  in  reasonable  agreement 
with  the  S  wave  section,  though  the  predicted  times  for  Sn  may  be  a  little  early. 
The  general  character  of  the  P  and  S  fields  is  similar  but  in  this  case  the  ratio  of 
the  mantle  to  crustal  amplitudes  is  larger  for  S  than  P. 

In  order  to  see  whetner  the  pattern  of  propagation  of  the  mantle  phases  is 
consistent  from  region  to  region  and  also  to  investigate  the  evolution  of  the 
Sn/Pn  amplitude  ratio  with  distance  we  have  searched  for  refraction  recordings 
with  good  S  wave  arrivals  in  the  distance  range  from  100  -  1000  km.  Although 
much  of  the  long-range  profiles  in  Europe  were  originally  recorded  with  three- 
component  sensors,  most  of  the  interpretations  have  been  based  on  the  P  wave 
data  from  the  vertical  component  records  and  relatively  few  S  wave  sections 
have  been  published.  However,  sufficient  S  wave  data  exists  for  us  to  be  able  to 
begin  to  make  a  comparison  of  the  propagation  characteristics  of  the  P  and  S 
wave  field. 

As  an  example  we  present  in  figure  6  a  comparison  of  P  and  S  recordings  at 
ranges  from  200  -  800  km  (Him,  1977)  from  a  profile  of  instruments  deployed 
across  France  to  record  a  large  explosion  fired  in  the  sea  off  Brest.  The  record 
sections  are  phased  to  concentrate  attention  on  the  mantle  arrivals  which  in 
nuclear  monitoring  work  we  would  normally  term  Pn  and  Sn.  However  we  can 
see  clearly  from  figure  6  that  such  phases  have  significant  substructure  and  the 
almost  constant  apparent  velocity  is  built  up  from  a  sequence  of  en-echelon 
phases  which  progressively  become  first  arrivals  as  the  distance  increases.  The 
interpretation  of  such  arrivals  in  terms  of  a  one-dimensional  velocity  model 
leads  to  very  complex  velocity  distributions  with  multiple  low-velocity  zones 
(see  e.g.  Kind  1974),  but  it  is  very  likely  that  much  of  the  complexity  arises  from 
lateral  heterogeneity  in  the  structure  in  the  upper  most  mantle  (Fuchs  &  Schulz, 
1976).  A  similar  pattern  of  complex  structure  within  the  wavegroup  characterised 
by  Pn  can  be  seen  in  the  record  sections  presented  by  Yegorkin  &  Chernyshov 
(1983)  from  peaceful  nuclear  explosions  in  the  Soviet  Union  (fig  13). 

In  figure  7  we  illustrate  the  pattern  of  rays  associated  with  the  class  of  model 
which  has  been  proposed  to  explain  the  complexity  of  the  observed  mantle 
arrivals.  A  weak  gradient  zone  below  the  Moho  is  terminated  by  a  low  velocity 
zone  (which  is  needed  to  get  a  suitable  delay  for  arrivals  from  greater  depth). 
This  leads  to  a  relatively  weak  Pn  arrival  from  the  uppermost  mantle  which  is 
supplan  .ed  at  greater  distance  by  the  P|  phase  returned  from  depth,  and  at  even 
greater  langes  by  Pn.  The  observations  of  Gajewski  et  al  (1990)  relate  to  the  Pn,  Sn 
arrivals  from  the  uppermost  mantle  and  not  to  these  arrivals  from  greater  depth. 
The  arrivals  Pj,  Pn,  Si,  Sn  are  frequently  more  energetic  and  are  the  phases  which 
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Figure  5:  P  and  S  wave  sections  for  the  vertical  component  of  displacement  for  an 
earthquake  in  northwestern  Scotland  recorded  during  the  long-range  refraction 
profile  through  the  British  Isles.  The  reduction  velocity  for  the  P  section  is  8 
km/ s,  and  for  the  S  section  it  is  4.62  km/s  (after  Kaminiski  et  al  1974). 
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Figure  6:  Three  component  record  sections  for  a  large  explosive  shot  off  the  French 
coast  recorded  across  a  long  range  profile  in  France  (after  Him  1977).  Reduction 
velocity  for  the  Z  section  is  8  km/s,  whereas  for  the  L  and  T  sections  the 
velocity  is  4.62  km/s 
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would  be  designated  in  discrimination  work  as  Pn,  Sn  because  an  individual 
seismogram  is  not  sufficient  to  recognise  the  substructure  of  the  wavefield. 

In  general  the  quality  of  S  data  from  explosive  sources  is  not  very  good  because 
the  P  coda  is  energetic  and  so  we  are  looking  for  S  arrivals  against  a  relatively 
noisy  background.  However  some  experiments,  especially  in  Precambrian 
terrains  where  the  shear  velocities  are  relatively  high,  do  give  a  very  good 
representation  of  the  S  wavefield  even  on  vertical  component  instruments. 
Figure  8  shows  the  locations  of  three  refraction  profiles  in  Fennoscandia  which 
have  good  quality  S  arrivals. 

Figure  9  shows  P  and  S  wave  sections  from  the  same  shotpoint  point  A  in  a 
Finnish-Polish  refraction  experiment  SVEKA  across  part  of  the  Baltic  Shield  in 
southern  Finland  (Luosto  et  al  1984,  Grad  &  Luosto  1987).  The  S  arrivals  are  very 
clear  and  give  a  good  general  agreement  with  the  pattern  of  the  P  arrivals.  In  this 
case  the  crust  is  rather  thick  and  the  transition  from  Moho  reflections  to 
refractions  occurs  at  the  end  of  the  profile.  Interestingly  the  reflection  from  the 
Moho  around  160  km  offset  is  clearer  on  the  S  wave  section.  Rg  is  only  visible  to 
100  km. 

The  quality  of  the  S  records  in  figure  9  is  surpassed  by  the  remarkable  set  of 
records  from  a  further  Finnish-Russian  experiment  BALTIC  which  lay  further  to 
the  southeast  (closer  to  the  Russian  border).  The  data  from  shot  point  B  is  shown 
in  figure  10  and  at  first  sight  is  difficult  to  distinguish  which  is  the  P  and  which  is 
the  S  profile  (other  than  by  examining  the  timing  of  the  Rg  phase).  The 
correspondence  of  the  phases  is  exceptionally  dose  and  here,  at  least,  there  is  a 
very  distinct  Sn  arrival  beyond  200  km.  The  reversed  profile  from  shotpoint  G 
has  less  clear  Pn  and  Sn  arrivals  with  a  cross  over  around  3020  km.  Despite  the 
pronounced  topography  on  the  Moho  (see  e.g.  Bannister,  Ruud  &  Husebye  1991), 
good  mantle  arrivals  are  seen  on  these  profiles. 

A  further  experiment  in  northern  Finland  and  Sweden,  the  POLAR  profile 
(Luosto  et  al  1989)  shows  uniformly  good  S  propagation  from  all  shot  points 
including  those  with  fairly  small  charges.  A  record  section  from  shotpoint  A  at 
the  southern  end  of  the  line  is  shown  in  figure  11.  In  this  area  it  would  appear 
that  the  Moho  is  very  sharp  because  of  the  very  distinct  PmP.  SmS  arrivals  which 
can  be  seen  to  connect  directly  into  the  Pn,  Sn  phases  at  greater  distance.  This 
profile  lies  relatively  close  to  the  Arcess  array  for  which  there  are  generally  very 
clear  observations  of  S  phases  from  the  mine  blasts  on  the  Kola  peninsula.  This 
geologic  province  is  the  oldest  in  Fennoscandia  and  it  wouiu  appear  that 
attenuation  is  rather  low. 

A  comparison  of  the  three  sets  of  P  and  S  record  sections  presented  in  figs  9-1 1 
shows  noticeable  differences  in  the  character  of  the  wavefield.  For  the  SVEKA 
profile  (fig  9)  there  are  strong  direct  Pg  and  Sg  phases  extending  out  to  200  km, 
where  they  are  overtaken  by  the  Pn  and  Sn  phases.  The  reflection  from  the  crust 
mantle  boundary  is  more  distinct  for  S  than  for  P,  indeed  it  is  quite  difficult  to 
correlate  the  P  wave  arrival  to  distances  smaller  than  180  km  where  the 
corresponding  S  reflection  is  very  clear. 

For  tl'e  BALTIC  profile  (fig  10)  which  is  again  in  southern  Finland,  Pg  and  Sg 
are  cleai  but  by  comparison  less  energetic  than  for  SVEKA.  In  this  case  both  the 
prograd  j  and  retrograde  travel  time  branches  associated  with  reflection  and 
refraction  at  the  crust-mantle  boundary  can  be  clearly  seen  for  both  P  and  S.  The 
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Fig.  1.  DSS  profiles:  6  -  Finlap;  c  -  Barents  Sea;  d- Pechenga-Umboiero;  t  -  Pechenga-Kostamuksha;  /-  Kem-Tulos;  g- 
Baltic;  A-Sveka;  /  -  Sylen-Porvoo;  J~  Ladoga;  A -Polar;  /- Kem-Uchta;  m  -  Elijarvi-Lahnaslampi.  SF- Svecofennian  Do¬ 
main;  KA  -  Karelian;  fl£  -  Belomordian  Province;  KO  -  Kola  Peninsula  Province;  CA  —  Calcdonides;  LGB  -  Lapland  Cranulite 
Belt;  LBB  -  Ladoga- [lothnian  Bay  Zone;  W-Viborg  Rapakivi  Granite  Intrusion  (tectonogeological  map  by  Korsakova  et  al., 

1987). 

Figure  8:  Configuration  of  the  SVEKA  and  BALTIC  and  POLAR  refraction  experiments 
in  Fennoscandia. 


figure  9:  Comparison  of  P  and  S  record  sections  from  the  SVEKA  refraction  profile  in 
Finland  using  vertical  component  seismometers,  the  traces  are  amplitude 
normalised.  The  reduction  velocity  for  the  P  section  is  8  km/s,  and  for  the 
S  section  it  is  4.62  km/s  (after  Luosto  et  al  1984,  Grad&  Luosto  1987). 


lorthem  Fennoscandia  using  vertical  component  seismometers,  the  traces 
amplitude  normalised.  The  reduction  velocity  for  the  P  section  is  8  km/s,, 
for  the  S  section  4.62  km/s  (after  Luosto  et  al  1989). 


strength  of  the  mantle  arrivals  near  the  crossover  at  180  km  leads  to  strong 
suppression  of  the  Pg  and  Sg  arrivals  in  these  amplitude  normalised  displays. 

For  the  POLAR  profile  (fig  11)  the  exceptional  strength  of  the  reflected  energy 
from  the  Moho  suppresses  the  display  of  Sg,  in  particular  between  120  and 
160  km.  If  a  phase  association  were  to  be  made  in  this  case  on  the  maximum 
S  amplitude  it  would  be  displaced  by  2-3  seconds  from  the  assumed  trajectory  for 
lg  and  could  well  lead  to  significant  mislocation  of  an  event.  The  mantle  arrivals 
are  weaker  than  for  the  BALTIC  profile  and  we  can  see  the  Pg  and  Sg  arrivals  grade 
into  the  PmP  and  SmS  arrivals  as  would  be  predicted  by  standard  crustal  models. 

The  map  of  refraction  profiles  in  Fennoscandia  indicates  a  number  of  profiles 
within  the  Soviet  Union,  with  a  relatively  dense  shotpoint  spacing. 
Unfortunately  almost  all  the  Deep  Seismic  Sounding  experiments  carried  out 
within  the  USSR  have  used  analogue  recordings.  Relative  close  station  spacing 
e.g.  100  m  with  a  5  km  recording  spread  has  enabled  close  correlation  of  features 
on  the  seismograms.  Attention  has  largely  been  focussed  on  the  time  relations  of 
the  different  features  and  in  published  accounts  only  very  limited  descriptions  of 
the  amplitude  behaviour  are  provided.  It  is  therefore  quite  difficult  to  extract 
systematic  information  on  the  character  of  the  wavefield. 

A  very  interesting  experiment  was  carried  out  by  Jentsch  (1979).  He  hand 
digitised  a  very  large  number  of  Russian  recordings  from  a  DSS  profile  in  the 
southern  Ukranian  shield,  and  then  attempted  to  apply  conventional  one¬ 
dimensional  interpretation  techniques  to  this  data.  Figure  12  displays  a  portion 
of  the  record  section  which  displays  both  the  P  and  S  energy.  The  source  for  the 
recordings  were  multiple  explosive  charges  in  separate  holes.  A  notable  feature 
of  these  records  is  the  relatively  short  correlation  distance  for  the  amplitude 
along  individual  phases.  Although  general  trends  are  clear,  there  are  significant 
variations  in  amplitudes  over  distances  little  larger  than  the  aperture  of  a 
regional  array  of  the  Noress/Arcess  type;  these  often  correlate  between  the  P  and 
S  phases  but  do  not  always  coincide  in  range.  We  therefore  certainly  need  to  be 
cautious  in  interpreting  amplitude  ratios  between  different  phases  from 
individual  seismograms.  For  discrimination  purposes  the  stability  of  an 
amplitude  ratio  estimate  based  on  a  number  of  closely  spaced  stations  is  certainly 
to  be  preferred  to  even  a  single  three  component  station.  Mini-arrays  with  a  few 
vertical  component  sensors  surrounding  a  central  three-component  systems 
would  seem  highly  desirable  in  this  context. 

Yegorkin  &  Chernyshov  (1983)  have  presented  record  sections  for  P  waves 
along  DSS  profiles  in  northern  Siberia  (to  1700  km)  and  central  Asia  (to  900  km). 
The  typical  spacing  of  traces  is  to  10  km  and  a  portion  of  the  Siberian  line  is 
reproduced  in  fig  13.  The  patterns  of  P  wave  propagation  for  the  onset  of  the 
seismogram  are  similar  to  those  seen  on  long-rangf  refraction  profiles  in  Europe. 
The  energy  return  from  the  mantle  consist  of  a  number  of  distinct  arrivals  which 
can  be  correlated  over  substantial  distances.  It  would  appear  that  the  recording 
interval  used  in  these  very  long  range  refraction  experiments  was  too  short  to 
include  the  S  arrivals  from  the  mantle  as  well. 

THEORETICAL  STUDIES 

In  association  with  the  observational  studies  we  have  also  considered  ways  in 

which  the  behaviour  of  the  Pn  and  Sn  can  be  modelled  with  allowances  for  the 
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effect  of  lateral  heterogeneity.  As  we  have  seen  above,  we  cannot  neglect  the 
influence  of  free  surface  multiples  and  so  we  cannot  easily  separate  the 
modelling  of  mantle  phases  from  the  problem  of  wave  propagation  in  the  crustal 
waveguide. 

For  Sn  we  can  adapt  the  coupled  mode  scheme  previously  used  by  Kennett 
(1986)  for  modelling  Lg  wave  propagation  in  the  crust  with  allowance  for 
coupling  to  Sn,  by  simply  enlarging  the  number  of  modes  considered  and  thus 
the  phase  velocity  span  which  can  be  modelled.  For  propagation  up  to  3-4  Hz  this 
would  require  coupling  calculations  spanning  up  to  200  modes  (i.e.  200  coupled 
complex  differential  equations).  Such  calculations  are  feasible  for  many 
frequencies  on  a  supercomputer  but  are  rather  time  consuming,  and  even  so  we 
are  restricted  to  slight  perturbations  of  the  crust-mantle  boundary. 

In  principle,  the  mode  coupling  procedure  can  also  be  extended  to  the  Pn  and 
Pg  parts  of  the  wavefield  but  this  would  require  more  than  a  1000  modes  for 
frequencies  higher  than  2  Hz  and  as  such  does  not  present  an  attractive  route.  If 
full  account  is  taken  of  the  possibility  of  reflection  we  would  need  to  solve  a 
differential  equation  for  a  1000x1000  complex  matrix:  even  with  the  most 
efficient  algorithms  this  would  be  very  difficult  on  current  computers. 

We  have  therefore  been  looking  at  alternative  representations  of  the  seismic 
wavefield  which  exploit  the  idea  of  coupling  between  different  horizontal 
wavenumbers  (which  underlies  the  mode  coupling  procedure)  but  which  use  an 
alternative  expansion  of  the  field  in  orthogonal  functions.  A  suitable  basis  set 
exists  if  the  free  surface  can  be  neglected  and  so  could  be  used  for  propagation  to 
around  300  km.  However  no  satisfactory  solution  has  yet  been  found  to  ailow 
inclusion  of  surface  reflections.  We  will  continue  to  investigate  methods  for 
modelling  the  propagation  of  the  mantle  phases  over  the  next  year. 

One  aspect  of  the  character  of  Pn  and  Sn  which  has  proved  amenable  to 
modelling  is  the  influence  of  topography  at  the  crust-mantle  boundary.  In  the 
Finnish  example  above  there  is  substantial  topography  at  the  Moho,  and  the 
long-range  refraction  profile  along  the  length  of  the  British  Isles  (Bamford  et  al, 
1976)  has  also  shown  significant  topography  at  this  level.  For  high  frequency 
waves,  an  approximation  to  the  effects  of  an  irregular  boundary  can  be  made  by 
introducing  a  generalised  transmission  coefficient  whereby  a  single  incident 
plane  wave  couples  into  a  spray  of  plane  waves  on  the  other  side  of  the  interface. 

This  approach  can  be  introduced  into  a  modification  of  the  reflectivity  method, 
in  which  ray  techniques  are  used  to  transfer  energy  from  the  source  to  the 
irregular  boundary  and  from  that  boundary  to  the  receiver;  but  a  full  calculation 
is  made  in  a  horizontally  layered  structure  lying  beneath  the  crust-mantle 
boundary.  However,  the  amplitude  kernel  in  the  plane  wave  response  for  this 
region  is  modulated  by  the  redistribution  associated  with  the  generalised 
transmission  coefficients  on  entry  and  exist  from  the  mantle.  The  method 
assumes  that  the  top  of  the  reflection  zone  averages  to  a  flat  level  assumed  in  the 
reference  model;  this  can  give  rise  to  slight  inaccuracies  in  the  amplitude  and 
timing  of  waves  reflected  at  the  underside  of  the  interface.  The  dominant 
arrivals  corresponding  to  diving  waves  will  be  correctly  represented. 

The  emplitude  distance  behaviour  for  a  model  derived  from  the  Moho 
topography  of  Bamford  et  a!  (1976)  is  displayed  in  fig  14,  together  with  the  form  of 
the  crust-mantle  interface  for  propagation  from  north  to  south  and  from  south  to 
north.  These  amplitudes  are  compared  with  those  calculated  for  a  laterally 
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uniform  reference  model  and  have  been  derived  from  the  synthetic  seismogram 
sections  in  Figure  15,  which  have  been  calculated  using  a  phase  velocity  window 
from  8.1  to  10  0  km/s  in  all  cases.  This  window  corresponds  to  purely 
propagating  waves  in  the  layer  immediately  below  the  crust-mantle  interface  and 
allows  direct  comparison  between  flat  and  perturbed  topography  at  the  Moho. 

The  results  for  the  reference  model  (with  horizontal  layering  and  a  Moho 
depth  of  32.5  km),  are  indicated  in  fig  14  by  the  solid  circles,  open  squares  are 
used  for  propagation  from  north  to  south,  and  open  triangles  are  used  for 
propagation  from  south  to  north.  The  fluctuations  introduced  by  the  presence  of 
topography  are  significant  but  in  general  less  than  the  typical  variability  in 
observed  amplitudes.  The  offset  distance  from  the  Moho  to  the  recording  point 
at  the  surface  is  about  45  km,  so  the  amplitude  symbol  will  be  displaced  from  the 
relevant  topography.  The  major  change  from  the  reference  model  is  a  shift  of  the 
amplitude  peak  near  300  km  to  shorter  distances  and  an  enhancement  of  the 
amplitude  minimum  near  400  km  range. 

These  features  are  clearly  seen  in  the  record  sections  of  the  theoretical 
seismograms  in  fig  15.  The  effect  of  topography  is  to  give  a  pronounced 
concentration  of  energy  near  300  km  accompanies  by  reduced  amplitudes  near 
400  km.  If  we  were  to  interpret  the  NS  and  SN  record  sections  in  terms  of 
horizontally  stratified  models,  we  can  see  that  the  resulting  velocity  models 
would  differ  significantly  from  the  reference  models.  In  particular,  to  change  the 
peak  in  Pi  amplitudes  a  more  pronounced  velocity  inversion  would  be  required. 

These  calculations  show  the  effect  of  topography  on  the  P  wavefield,  and  we 
can  anticipate  similar  effects  for  S  waves.  However,  unless  the  velocity  ratio 
between  P  and  S  waves  is  a  constant  throughout  the  medium,  the  propagation 
paths  for  P  and  S  waves  will  differ.  The  focussing  and  defocussing  effects 
associated  with  topography  can  therefore  be  displaced  slightly  between  the  P  and 
S  wave  patterns  which  will  give  rise  to  fluctuations  in  the  ratios  of  amplitudes  of 
mantle  P  and  S  waves  and  so  complicate  use  as  a  discriminant. 

CONCLUSIONS 

We  have  been  able  to  demonstrate  that  the  amplitude  behaviour  of  Pn  and  Sn 
with  distance  is  quite  complex  and  that  it  is  necessary  to  look  at  the  anatomy  of 
the  phases  in  order  to  understand  the  potential  performance  of  some  regional 
discriminants.  In  detail  the  Pn  and  Sn  phases  tend  to  break  up  into  a  sequence  of 
overlapping  branches  associated  with  fine  structure  in  the  upper  part  of  the 
lithosphere.  Anomalously  low  Sn  amplitudes  can  arise  between  200  and  300  km 
away  from  the  source. 

The  current  class  of  regional  phase  discriminants  based  on  averaged  behaviour 
of  seismic  phases  such  as  Pn  and  Sn  could  give  misleading  results  in  some 
distance  ranges  because  of  the  way  in  which  the  peak  amplitudes  in  the  requisite 
time  windows  shift  between  different  wave  groups.  Such  effects  are  difficult  to 
pick  up  from  isolated  earthquakes  or  mine  blasts  but  are  evident  in  long  range 
refraction  profiles  or  comparable  rec  ord  sections. 

It  will  be  difficult  to  assess  the  amplitude  behaviour  of  the  Pn  and  Sn  phase  for 
the  area  around  a  particular  receiver  site  without  careful  calibration,  and  the 
analysis  of  regional  phases  will  need  to  be  adapted  to  take  the  nature  of  the 
particular  phases  into  consideration.  The  maximum  amplitude  within  a 

21 


specified  group  velocity  window  may  not  be  the  best  measure  of  the  P  or  S  wave 
content  of  the  seismogram. 
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